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Abstract. We have studied quasielastic charged current hyperon production induced by on free nucleon and the nucleons 
bound inside the nucleus and the results are presented for several nuclear targets like 40 A r, 56 Fe and 2m Pb. The hyperon- 
nucleon transition form factors are determined from neutrino-nucleon scattering and semileptonic decays of neutron and 
hyperons using SU(3) symmetry. The nuclear medium effects(NME) due to Fermi motion and final state interaction(FSI) 
effect due to hyperon-nucleon scattering have been taken into account. Also we have studied two pion production at threshold 
induced by neutrinos off nucleon targets. The contribution of nucleon, pion, and contact terms are calculated using Lagrangian 
given by nonlinear a model. The contribution of the Roper resonance has also been taken into account. The numerical results 
for the cross sections are presented and compared with the experimental results from ANL and BNL. 
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INTRODUCTION 

The neutrino oscillation experiments being done using accelerator (anti)neutrino beams in the energy region of few 
GeV are also providing cross section measurements of various reactions induced by neutrinos and antineutrinos on 
nuclear targets which are needed for validating various Monte Carlo neutrino event generators. The reported cross 
sections are available mainly for elastic, quasielastic and single pion production processes on 12 C and 16 O nuclei. 
The availability of high intensity (anti)neutrino beams in present generation neutrino oscillation experiments has 
opened up the possibility of cross section measurements of other processes like production of hyperons, strange 
baryons and multiple pions in the energy region of a few GeV in the experiments being done at MINERvA [1], 
MINOS [2], NOvA [3], LBNE [4], T2K [5] and MiniBooNE [6], In the lower energy region relevant for T2K [5] and 
MicroBooNE [7] experiments, any observation of these processes will help to understand the reaction dynamics of 
the threshold production of strange mesons and two pions. In view of this scenario, some calculations have been done 
recently for the production of hyperons [8], kaons [9, 10] and two pions [11], 

In this work we present a brief overview of the quasieleastic production of hyperons induced by antineutrinos on 
nucleons and nuclei and charged current two pion production induced by neutrino from nucleon. First, we are going to 
present the formalism, results and discussion for the quasieleastic production of hyperons from nucleons and nuclear 
targets followed by the formalism of two pion production from nucleons and their results and discussion. 


QUASIELASTIC PRODUCTION OF HYPERONS 
Formalism 

The quasieleastic production of hyperons induced by antineutrinos has been studied in the past using Cabibbo 
theory with SU(3) symmetry [12] as well as in quark models [13]. In addition to providing neutrino-nucleus cross 
section to be used in modeling the MC neutrino event generators, the differential cross sections also provide an 
opportunity to study nucleon-hyperon transition form factors at high Q 2 , which are available only at low Q 2 through 


the analysis of Semileptonic Hyperon Decays(SHD) [14]. An independent determination of these form factors help to 
test various assumptions made in the analysis of these processes like SU(3) symmetry, G invariance, absence of Flavor 
Changing Neutral Currents(FCNC), A Q = AS rule. Conserved Vector Current(CVC) and Partial Conservation of Axial 
Current(PCAC) hypothesis in the Standard Model(SM) when it is extended to the strangeness sector. Since these 
hyperons decay primarily through pionic decay modes, they also contribute to the pion production processes induced 
by antineutrinos, which is dominated by the pion production through A-excitation. In the low energy region where 
A-excitation is inhibited by threshold effects, the Cabibbo suppressed pion production through hyperons (A and E) 
production may be important. 

We consider the following processes 


Vfi(k)+p{p) p + {k')+A{p') 

Mk)+p(p) ->/i + (*0+zV) 

Vju(k) +n(p) -> p + (kf)+I.-(p'), (1) 

where k(k’) and p(p') are the momenta of intialffinal) lepton and nucleon. The differential scattering cross section 
is given by. 
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The form factors fj(q 2 ) and gi(q 2 ) are determined using T invariance, G invariance, SU(3) symmetry (symmetry 
properties of weak currents) like CVC and PCAC hypothesis. These symmetry considerations yield [12]: 

1. fi(q 2 ) and 8i(q 2 ) as real quantities 

2- h{q 2 ) = gi{q 2 ) =0 
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5. The form factors, f\p(q 2 ) = ^(//^(tf 2 ) +2/" 2 (? 2 )) and f\,i{q 2 ) = ~\Jif\2^q 2 ) respectively for transitions 
p —¥ E° and p — > A. Similarly gi(q 2 ) is -^%-ppgA{q 2 ) f° r P — ► E° and is s A (q 2 ) for p — » A transitions. 

We have used D = 0.804 and F = 0.463 and the following q 2 dependence of the electroproduction and weak form 
factors [15]: 
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The Sach’s form factors G E n (q 2 ) and G p ^ l {q 2 ) are parameterized as 
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The numerical values of various parameters are taken as 
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For gAifl 2 ) a dipole form has been taken i.e. gA(q 2 ) 
mass Ma = lGeV. 


with ‘ ?A ^ 


1.267 and the axial dipole 


Nuclear Medium Effects and Final State Interaction Effects 


When these reactions take place on bound nucleons in nuclear medium, Fermi motion and Pauli Blocking effects 
of nucleons are to be considered. In the final state after hyperons are produced, they may undergo strong interaction 
scattering processes through charge exchange) E~p —» An, Ap — > E + n, etc.) and inelastic (AN T,°N, YN —> Y'N') 
reactions like changing the relative yield of E°, E and A produced in the initial reactions shown in Eq. 1. In a special 
case, E + will appear as a result of final state interaction which are initially not produced through v p N —> jl + N reaction 
due to A Q = AS rule. 

The nuclear medium effects are calculated in a relativistic Fermi Gas model using local density approximation and 
the nuclear cross section is written as 
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where n^(p, r) is local occupation number of the initial nucleon of momentum p and is 1 for p < pf N and 0 otherwise 
with 

Pf„ = [3?rp„(r)] 1/3 for N = n,p (12) 

The final state interaction of hyperon-nucleon system is calculated in a Monte Carlo simulation approach. In this 
approach an initial hyperon produced at a position r within the nucleus which interacts with a nucleon to produce a 
new hyperon-nucleon state / = YfNf, within a distance / with probability P Y dl where P Y is the probability per unit 
length given by 

Py = a Y+n -yf{E) p n {r) + a Y+P -rf{E) Pp{r), (13) 

where p„(r)(p p (r)) is the local density of neutron (proton) in the nucleus and a is the cross section for YN —> f 
process. Now a particular channel / is selected. For the selected channel /, a hyepron-nucleon state is chosen by 
randomly selecting the momentum of initial nucleon consistent with Pauli blocking. A random scattering angle is 
generated using isotropic cross section for the hyperon-nucleon scattering cross section a, which then determines 
momentum(energy) of final hyperon and nucleon. If the momentum of final nucleon is above Fermi momentum, a 
new final state of YN system is obtained. This process is continued until the hyperon gets out of the nucleus. All the 
channels i.e. Y + n —> Yr + Nf and Y + p —► Yp + Np leading to a final state / are considered. These hyperons decay 
into pions through the processes. 

E°(A°) pn^pin 0 , E + pn°,nn + , E - nn~ (14) 

In the final state, a p 1 will be accompanied either by a K° or n and rarely by a K which will be produced as a 
result of the final state interaction in the nucleus. Note that it can also be produced as final state interaction effect when 
K°p —> K + n reaction takes place. 


Results 

In Fig. 1, we have presented the results for the total cross section for V p + p —> p - A reaction for free nucleons 
and compared them with the experimental results [16, 17, 18, 19, 20, 21]. We have also compared our results with 
the other theoretical results in quark model obtained by Wu et. al. [13] and Finjord et. al. [22] and in Cabibbo model 




FIGURE 1. c vs Ey , for Vfi + p — > /j + + A process. Experimental results (triangle right [16], square [17], triangle up [18], 

circle [19]), triangle down(cr = 2.6^2'j x 10 _40 c/« 2 ) [20] are shown with error bars. Theoretical curves are of Erriquez et 
al. [17](dashed-double dotted line), Brunner et al. [19](dashed line), and Kuzmin and Naumov [21](double dashed-dotted line) 
obtained using Cabibbo theory with axial vector dipole mass as 1 GeV, 1.1 GeV and 0.999GeV respectively, while the results of 
Wu et al. [13](dotted line) and Finjord and Ravndal [22](dashed dotted line) are obtained using quark model. The results of present 
calculation are shown with solid line. Notice that we have multiplied the results of Wu et al. [13] by 3 to plot on the same scale. 


by Kuzmin et. al. [21]. We see that the most recent results in quark model by Wu et. al. [13] underestimate the cross 
sections. The theoretical results presented in Fig. 1 do not include nuclear medium effects while the experimental 
results are on nuclear targets. Obviously experimental data with better statistics are needed to study nuclear medium 
effects. However, Cabibbo theory with SU (3) symmetry seems to work well for explaining the present results on 
Vjjp — > Li A. In Fig. 2, we show the results for v^p —> p L° where new data are needed to draw any conclusion. For 
completeness we present the results for v^p —> id 1 A in 40 ,4r. 5( Fe and lm Ph in Fig. 3 where we also show the effect 
of final state interactions. 



FIGURE 2. a vs Ey . for v^+p —> p + +E° process. Experimental points is taken from [17]. Present results are shown with 
solid line. Also the results of Wu et al. [13](dotted line) and Finjord and Ravndal [22](dashed dotted line) have been presented. 
Notice that we have multiplied the results of Wu et al. [13] by 3 to plot on the same scale. 


In Fig. 4, we present the results for 

f or 56 p„ 208 

spectrum i.e. 


do 

W- 


in 40 Ar averaged over the MiniBooNE antineutrino spectrum [6] and 


Fe and Z K Pb averaged over the MINERvA antineutrino spectrum [23]. It should be noted that lepton energy 
^ can be easily obtained from the Q 2 distribution in each case. The details are given in Ref. [15]. 










FIGURE 3. a vs E ^ for 40 Ar, 56 Fe and 2m Pb nuclei. 



FIGURE 4. < > vs 2 2 (Left panel) in 12 C, 56 Fe and 208 Pft nuclear targets obtained by averaging Q 2 -distribution over the 

MINERvA [23] flux for the reaction given in Eq. 1. The results are presented with nuclear medium and final state interaction 
effects. 


THRESHOLD PRODUCTION OF TWO PIONS 


Two pion production is the threshold process of inelastic reactions beyond single pion production. Experimental 
observation of following processes on nucleon targets 
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have been reported at ANL [24] and BNL [25] experiments. To the best of our knowledge no theoretical calculation 
has been reported in literature except in the threshold region [11, 26]. These calculations make use of an effective 
chiral Lagrangian to describe the interaction of weak currents with pions and nucleons. It seems that the resonance 
contribution is as important as the contribution of nucleons and pions(including the contact term) even in the threshold 
region. It may turn out that the 2n production is dominated by processes in which a pion is produced along with a A in 
the intermediate state giving rise to another pion as happens in the electroproduction [27, 28]. Calculations of two pion 
production in an intermediate A dominance model will be very useful in analyzing the existing data from ANL [24] 
and BNL [29] experiments on two pion production and data to be obtained in future from the experiments looking for 
neutrino oscillation. 



















































FIGURE 5. Nucleon pole, pion pole and contact terms contributing to 2k production. 


In this work, we report the calculation of weak charged current production of two pions in threshold region. We 
use a nonlinear chiral Lagrangian model which has been extensively used to study the single pion production [30], 
The Lagrangian describes the interaction of pions and nucleons and generates the week vector and axial vector current 
under a SU{ 2) x SU( 2) transformation which interacts with W ± IZ boson as described by the Standard Model. The 
effective Lagrangian is given as 
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where ^ ^ J is the nucleon field, (f is the isovector pion field, T are the Pauli matrices and f K = 93 MeV 

is the pion decay constant. The vector(Vu) and axial vector)//,) currents generated by this Lagrangian under chiral 
transformation are given by 
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These currents couple to VV' + (Z) boson for charged current) neutral current) interaction with nucleon and pion as 
described in the Standard Model. Various terms in Vu and A ^ describe the vector and axial vector couplings for 
WNN, WNn WNNn etc. vertices. And the matrix elements for various Feynman diagrams can be calculated for 
vN —> j,L nnN process using the rules of covariant perturbation theory. Using these currents given in Eqs. 19 and 20 
we obtain following 16 diagrams for 2n productions shown in Fig. 5, for which the matrix element can be explicitly 
written using W + +N —> N, W + + N —> Nn, W + + N —> Nnn and W + n —> n vertices from Eqs. 19 and 20. The 7ZNN 
and kkNN etc. couplings are derived from Eq. 18. With these matrix elements for the hadronic transition current , 
the matrix element for the process ( k ) + N(p) —> p~{k’) + n(k Xl ) + n(k Xl ) +N(p') is written as 
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where 

^ =u{k')f(\-r’)it{k) (22) 

and 

Jn = {■N*(p')n(k Kl )n(k* 2 )\jn\N(p )> (23) 

The present formalism gives the correct form of the nucleon vector (V?) and axial vector (A?) currents which 
couple to the W ± boson but go with a point coupling for the WNN vertex. The form factors are phenomenologically 
introduced at this vertex which are consistent with electron proton and neutrino nucleon scattering. The following 
forms are used for the matrix element for VT and A v currents 
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with fj 2 (<? 2 ) = f" 2 (q 2 ) where the vector form factors f^ 2 and the axial vector form factor gi {q 2 ) are given 

by Eq. 8. 

The explicit expressions for the matrix elements of 16 Feynman diagrams shown in Fig. 5 contributing to are 
given in Ref. [11]. The cross section is then expressed as 
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CONTRIBUTION OF iV*(1440) RESONANCE 

The Roper A* (1440) is the lowest lying resonance with significant coupling to two pion decay mode in S state 
and is expected to contribute in the threshold region. The N* —> Nnn coupling has been studied in nN —> nnN 





FIGURE 6. Direct (left) and crossed (right) Roper excitation contributions to 2 K production. 


and NN —> NNnn reactions [31] where N* contribution found to be important. The N* —> Nun is described by the 
Lagrangian 


,m~ _ -J 2 „ 1 
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with cf = -7.27 GeV" 1 , c| = 0 GeV" 1 [31]. 

The Feynman diagram for V^N —> u Nkk through N* excitation is shown in Fig. 6, for which the matrix element 
is written as 

JC = 2g*u{p')[SR{p + q)J%(q)+J%(q)S R {p'-q)\u{p) \ /£ = (29) 

where Sr is the propagator for Roper resonance 
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with /I = Af.v + Mr; M,y and Mr are respectively the mass of nucleon and Roper resonance and g* = —ct . 
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The axial vector form factor Gp is obtained using PCAC hypothesis as 
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with / is strength of N*Nn coupling determined by the N* Nn decay width F 
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with W as center of mass energy of nN state and |q cm | is the center of mass momentum. With F = 350 MeV, we obtain 
/ = 0.48. For G A {q 2 ), q 2 dependence is assumed to be of dipole form i.e. 
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with M,\ = 1 GeV. The isovector vector form factors F* = F* — F* n and F 2 = F 2 — F*„ are determined from the 
helicity amplitudes defined as 
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where the polarization vectors are given by 
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e ± = -^=(0,t1,-«,0), 

and for a photon of momentum q moving along the positive z-axis 
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with N for proton or neutron, a = 1/137, q is the momentum of the virtual photon, Icr = (W 2 — M 2 )/2W and W is the 
center of mass energy of the Roper. 

The EM y/V —> N* current is written as 
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Using Eqs. 36-39, we obtain the following relations: 

aN \/i = \q\g(<f) 


r F N *(n 2 

(< q a 4 - ? f) + i-- J —. —-<7“' q v 




irN* 

^2 




*(?)■ 


„2 jrN* 

q 


[ jU W+M jU- J 


> 7/2 = ^\^ 2 ) 


j^N* 

7i 1 " 


jrN* 

t 2 


i 


jtt W + M 


(39) 


(40) 

(41) 


with 


S(C) 


87T a(W+M)W 2 


y M(W-M)((W+M) 2 ~q 2 ) ’ 


(42) 


Inverting Eqs. 40 and 41 to obtain and f/ Vi in terms of ,4 / , and Sy 2 and using quark model predictions of 
A\ /2 = —2/3 Ay 2 and S"j 2 = 0 we obtain 
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The q 2 dependence of Ff* and F 2 * have been obtained by fitting the data on proton-Roper electromagnetic transition 
form factors to the experimental results in helicity amplitudes and their forms are obtained as 
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RESULTS 

We present numerical results for the different channels using the matrix elements corresponding to various Feynman 
diagrams shown in Fig. 5 and Fig. 6. All the coupling constants and form factors are fixed as described in the text 
through Eqs 18-20. The relative phases between different diagrams in Fig. 5 are fixed by the Lagrangian itself. 

In the case of resonance contribution the sign of N*Nn coupling is taken to be same as NNn coupling. The axial 
form factors are used as given in Eqs 31-33. For vector form factors F/ and F 2 we have considered four cases i.e. 
FF1, FF2, FF3, FF4 as follows. 

FF1: Our fit of helicity amplitude given by MAID [11] 

FF2: Form Factors as determined in the quark model of Meyer et al. [32] 

FF3: Parameters given by Lalakulich et al. [33] 




background and the N* contribution is not shown. See text for details. 


FF4: MAID analysis of helicity amplitudes [34] 

In Fig. 7, we present the results for the cross section for the process v u n —> j .1 pn ' n . We show separately the 
contribution of background terms(pion and nucleon, including the contact terms) and Roper resonance. We see that 
at lower energies the contribution of Roper is comparable to the background terms and can be more or less than 
background contribution depending upon the form factors used. The resonance contribution is sensitive to F% (q 2 )- 
At higher energies E v > 0.7 GeV, background terms dominate. In Fig. 8, we present the results for this process with 
experimental results [29, 35] in kinematical region restricted by 
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Obviously our model underestimates the experimental results indicating the contribution of other resonances or even 
the presence of intermediate An state in this kinematics region. 

In Fig. 9, we present the results for v^p -A jl nun' and compared them with the experimental results [35]. Here 
again, our results are lower than the experimental results in this channel. It should be noted that while our results 
underestimate the experimental results in all these channels, they are larger than the results obtained by Adjei et. 
al. [26]. In view of this an improved calculation is highly desired for weak pion production in threshold region as well 
as at higher neutrino energies. In Fig. 10 we present our prediction for the cross section for other channels like: 


VfiP -a 

jl^p7t + 7l° 

(49) 

V/iP -a 

ji~nn + n + 

(50) 

V* 

ji~nn + n° 

(51) 

V^n —>- 

ji~pn°n° 

(52) 

CONCLUSION 



In this work we have presented a review of the theoretical calculation of reactions for the production of leptons 
due to quasielastic production of hyperons v^N -A jl 1 Y and two pion production \\,N —> jl nn'N. These are the 
lowest threshold processes beyond the single pion production and may be seen in the present generation of neutrino 









FIGURE 8. Cross section for the v^n —> fl p7l + 7t with cuts as explained in the text. Dashed line: Background terms. Solid 
line: Full model with set 1 of nucleon-Roper transition FF. Data are from Ref. [29] (solid circles) and Ref. [35] (open squares). 



FIGURE 9. Cross section for the v^p —> p n7t + 7t + with cuts as explained in the text. Note that there are no contributions from 
the7V*(1440) resonance to this channel. Data are from Ref. [35]. 


oscillation experiments some of them which are producing results on cross section measurements. Phenomenological 
cross section have been used for these processes in modeling the neutrino cross section for validating the Monte 
Carlo generators for neutrino events. It is desirable that a theoretical calculation of the cross section for these inelastic 
processes is made available for modelings the neutrino-nucleus cross section in this energy region. 

A theoretical understanding of quasielastic hyperon production processes will help us to understand the role of 
symmetries of weak interaction currents at higher momentum transfer in the quasielastic production of hyperons. 
In the case of two pion production, there is an urgent need of going beyond the region of threshold production to 
understand the reaction dynamics of two pion production in weak interaction. Even in the threshold region, there is 
a need to include other contributions not considered in this work. Experimentally, efforts should be made to observe 
these reactions in the ongoing experiments looking for neutrino oscillation in the few GeV energy region. 






























E y (GeV) 


FIGURE 10. Cross sections as a function of the neutrino energy. All calculations correspond to the full model with the FF1 set 
of nucleon-Roper transition form factors. 
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